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Abstract: In this work we use the (/-oscillator formalism to construct the atypical (short) 
supersymmetric representations of the centrally extended U q (su(2\2)) algebra. We then 
determine the S-matrix describing the scattering of arbitrary bound states. The crucial 
ingredient in this derivation is the affine extension of the aforementioned algebra. 
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1 Introduction 

Integrable systems constitute a special class of models in mathematics and physics. Their 
properties allow them to be solved exactly and thus they form a very useful playground 
for studying various systems. One common feature shared by these models is that they 
are closely related to some underlying algebraic structures. Hence most of the quantum 
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integrable systems exhibit a large and powerful symmetry algebra, for example of Yangian 
or quantum affine type. A particularly interesting example is the Hubbard model. 

The Hubbard model, which was named after John Hubbard, is the simplest model of 
interacting particles on a lattice. It has only two terms in the Hamiltonian: the hopping 
term (kinetic energy) and the Coulomb potential [1]. The model describes an ensemble of 
particles in a periodic potential at sufficiently low temperatures such that all the particles 
may be considered to be in the lowest Bloch band. Moreover, any long-range interactions 
between the particles are considered to be weak enough and are consequently ignored. It 
is based on the tight-binding approximation of superconducting systems and the motion of 
electrons between the atoms of a crystalline solid. Despite its apparent simplicity, there are 
different applications and generalizations describing a plethora of interesting phenomena. 
In the case when interactions between particles on different sites of the lattice can not 
be neglected and are taken into account, the model is often referred to as the Extended 
Hubbard model. The particles can either be fermions, as in Hubbard's original work, or 
bosons, and the model is then referred as either the Bose-Hubbard model or the boson 
Hubbard model. The latter can be used to study systems such as bosonic atoms on an 
optical lattice (for a decent overview of various generalizations see reprint volumes [2-4] 
and also a more recent book [5]). 

A very specific class of models that share features with the one-dimensional Hubbard 
model and the supersymmetric t-J model [6] is the so-called Alcaraz and Bariev model 
[7]. It contains an extra spin-spin interaction term in the Hamiltonian and it shows some 
characteristics of superconductivity. This model can be viewed as a quantum deformation 
of the Hubbard model in much the same way as the Heisenberg XXZ model is a quantum 
deformation of the XXX model. This model has a specific R-matrix which can not be 
written as a function of the difference of two associated spectral parameters. This paradigm 
is related to the very interesting but at the same time complicated algebraic properties of 
the model. 

In recent years there has been renewed interest in integrable models arising from the 
discovery of integrable structures in the context of the AdS/CFT correspondence. For a 
recent review see [8] and references therein. The worldsheet S-matrix encountered there is 
one of the central objects of research and it turns out to have a lot in common with the 
specific cases of the Hubbard model considered in [9, 10]. Interestingly, the S-matrix of 
such a Hubbard model is obtained as a special limit of this worldsheet S-matrix [11]. 

The exact integrability of the one-dimensional Hubbard model was established by 
B. Shastry [12]. It was also shown that the model exhibits y(su(2)) © y(su(2)) Yangian 
symmetry [13]. However this symmetry is insufficient to constrain Shastry's S-matrix 
completely. Similarly, the worldsheet S-matrix for the AdSs x S 5 superstring also turns out 
to have Yangian symmetry [14] . However the Yangian in this model is based on a larger Lie 
algebra, the centrally extended su(2|2) Lie superalgebra. This underlying Lie superalgebra 
turns out to be powerful enough to constrain the S-matrix [15-17] (up to an overall phase, 
the so-called 'dressing factor' [18-20]) in the case where at least one of the representations 
is fundamental. However, Yangian symmetry (or equivalently the Yang-Baxter equation) is 
required in order to find the S-matrix describing the scattering of states that live in higher 
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representations [21, 22]. This specifically concerns the bound states in the system which 
transform in supersymmetric short representations [23-26]. The bound state scattering 
matrix can be explicitly constructed with the help of the underlying Yangian symmetry 
[27]. 

Nevertheless, there are still some problems concerning this infinite dimensional Yangian 
algebra due to some of its unusual features. The centrally extended su(2|2) Lie superalgebra 
has a degenerate Cartan matrix which prohibits the direct application of most of techniques 
related to the theory of Yangians. For the case at hand this has been partially circumvented 
in several ways: by enlarging the algebra by an sl(2) automorphisms [14], by considering the 
a — > limits of the exceptional Lie superalgebra 9(2, 1; a) [28] or building Drinfeld's second 
realization [29]. However this still proves to be an obstacle when, for example, one tries to 
construct the universal R-matrix [30, 31]. This object encodes all the scattering data in 
the theory in a purely algebraic form. Another issue that is not completely understood is 
the appearance of the so-called secret symmetry [32]. This is an additional symmetry of 
the S-matrix that does not have a corresponding Lie algebraic generator in the extended 
su(2|2) and could be interpreted as a u(l) outer automorphism of su(2|2). Resolving these 
issues could shed some light on the underlying algebraic structures and put the methods 
used to solve the model on a more firm footing. 

A possible route for attacking these issues was put forward in [10], where the quantum 
deformation Q of the extended su(2|2) algebra was studied. This g-deformed algebra has a 
number of interesting features such as a rather symmetric realization of the different central 
elements. Excitingly, just as in the non-deformed case, there is a link to Hubbard models, 
more specifically it describes spectrum of deformed supersymmetric one-dimensional Hub- 
bard models [10, 33]. The undeformed Hubbard model is revealed by taking a specific limit 
of deformed model [34]. Moreover, by sending the quantum deformation parameter q — > 1, 
the S-matrix under the consideration reduces to the AdS/CFT worldsheet S-matrix. As 
such, this matrix encompasses both different varieties of Hubbard models as well as the 
AdS/CFT worldsheet S-matrix and seems to provide a unifying algebraic framework for 
describing this class of models. 

The q-deformed S-matrix in the fundamental representation is constrained up to an 
overall phase by requiring invariance under Q itself. However, in the light that both 
the AdS/CFT and the Hubbard model S-matrices are actually invariant under an infi- 
nite dimensional symmetry algebra, it should not be surprising that such a structure is 
also present here. Indeed, the larger algebraic structure underlying this S-matrix is the 
quantum affine algebra Q [35]. This infinite dimensional algebra is obtained by adding an 
additional fermionic node to the Dynkin diagram of Q. In the q — >■ 1 limit one can retrieve 
the Yangian generators of centrally extended su(2|2) by considering the appropriate com- 
binations of generators of Q. This fuels the idea that Q plays a similar role as the Yangian 
in the undeformed case. More specifically, it is expected that the S-matrix in the higher 
representations is uniquely defined up to an overall phase by the underlying quantum affine 
algebra Q. This indeed turns out to be the case as we will show in this work. 

The class of representations we are considering in this work are the supersymmetric 
short representations. These representations are called short because the central elements 
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are not independent; they satisfy the so-called shortening condition [23]. In order to 
construct these representations, we employ the formalism of quantum oscillators. It is 
a quantum version of the well-known harmonic oscillator algebra and is defined as 



The use of quantum oscillators in the context of quantum groups was investigated earlier 
in [36-38]. By employing Fock space type modules, (/-oscillators naturally give rise to the 
representations of quantum groups. This approach was first formulated for the quantum 
deformed algebra U q (s[{2)) and later extended to simple Lie (super)algebras of a more 
general type, see e.g. [39]. Since then quantum oscillators have become an important part 
of the theory of quantum deformed algebras. 

Apart from being an interesting mathematical playground for studying the quantum 
affine algebra Q and its S-matrix, there is also a more elaborate motivation for considering 
these representations and the corresponding S-matrix. Firstly, there might be some possible 
applications in the context of the deformed Hubbard model. Secondly, it turns out that 
bound state states transform exactly in these representations of (/-oscillator algebra. It is 
important to study bound states for many reasons. For example, bound states usually play 
a crucial role in the thermodynamics of the model. In the case of the non-deformed model 
in AdS/CFT, the thermodynamic Bcthc ansatz (TBA) formalism is key in describing the 
complete spectrum of the theory [40-43]. The bound state S-matrix then governs the large 
volume solutions of both the TBA equations and the Y-system. Thus this is one of the 
first steps towards the TBA and Y-system formalism for the (/-deformed model. And, 
consequently, it might give some useful insights in these structures in the context of the 
AdS/CFT superstring. For example, there might be an interesting link to the recently 
constructed (/-deformed Pohlmeyer reduced version of the superstrings in the AdSs x S 5 
background [44, 45] which seems to be closely related to the (/-deformed model constructed 



In this work we derive the matrix structure of the general bound state S-matrix by 
employing the methods used in the context of the AdS/CFT superstring [27], but rather 
than using the Yangian symmetry we make use of the underlying quantum affine algebra 
Q. Our approach is based on the identification of invariant subspaces in the scattering 
theory that are specified by their invariance properties under the Cartan elements of the 
algebra. Then we use the rest of the algebra generators to relate these subspaces to each 
other resulting in the explicit form of the corresponding S-matrix. Just as in [27] we find 
the S-matrix in a factorized form reminiscent of the Drinfeld twist [46] . 

The paper is organized as follows. In section 2 we discuss the quantum deformation 
Q of the extended ZY(su(2|2)) algebra and its affine extension Q. Then in section 3 we 
introduce the quantum oscillator formalism and construct the supersymmetric short rep- 
resentations of Q. In section 4 we present the explicit derivation of the S-matrix for these 
representations. Subsequently, in section 5, we specify some explicit cases, we reproduce 
the fundamental R-matrix and also we give the precise form of the scattering matrix when 
one of the spaces forms a fundamental representation. We end with a brief discussion on 
the results and interesting directions for future research. The majority of the S-matrix 
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Figure 1. Dynkin diagram for the su(2|2) algebra. 

coefficients and results of the intermediate steps of the performed calculations are spelled 
out in the appendices. 

2 Quantum afRne algebra of extended U q (su(2\2)) 

In this section we review the quantum deformation of the extended su(2|2) algebra [10] 
and its affine extension [35]. 

2.1 Quantum deformation of extended su(2|2) 

The quantum deformed extended su(2|2) algebra Q was introduced in [10]. This algebra is 
generated by the three sets of Chevalley-Serre generators {Ej, Kj, Fj} (j = 1,2,3) where 
Ej and F~ are raising and lowering generators respectively and Kj = q H i are the Cartan 
generators. We will consider the case when E 2 and F 2 are fermionic generators and the 
rest are bosonic. This corresponds to the su(2|2) Dynkin diagram in Figure 1. In addition, 
this algebra has two central charges U and V = q and two parameters: the deformation 
parameter q and the coupling constant g. There is also a third parameter a, which describes 
the relative scaling of E2 and F2. Even though it is possible absorb this parameter into 
thes generators by a suitable redefinition, we will keep it unspecified. 

Algebra. The commutation relations which include the mixed Chevalley-Serre generators 
are (j, A; = 1,2,3) 

K 3 E k = q +DAjk EfrKj, KjF k = q~ DAjk F k Kj, [Ej,F k } = I),/,, ' <J ^\ , (2.1) 

where the associated Cartan matrix A and normalization matrix D are given by 

/ +2 -1 0\ 

DA =1—1 0+1, L> = diag(+1,-1,-1) . (2.2) 
\ 0+1-2/ 

There are also the unmixed commutation relations, called the Serre relations (j = 1,3), 

[E lt Ea] = {E 2 ,E 2 } = [Ej, [E V E 2 ]] -{q-2 + q~ 1 )E J E 2 E J = 0, 

[F u F 3 ] = {F 2 ,F 2 } = [Fj, [F V F 2 ]} - (q - 2 + q'^F^Fj = 0. (2.3) 

In addition, this algebra satisfies the extended Serre relations that give rise to two central 
elements U and V as follows, 

ga(l - U 2 V 2 ) = {[£ 2 ,M [£2, £3]} ~(q-2 + q-^E^E^ 
ga -\V- 2 -U- 2 ) = {[F 2 ,F 1 ],[F 2 ,F 3 ]} -( q -2 + q- 1 )F 2 F 1 F 3 F 2 . (2.4) 
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Figure 2. Dynkin diagram for the affine su(2|2) algebra. 

The central element V is also related to the Cartan generators through 

V- 2 = KiKlK z . (2.5) 
The conventional U q (su(2\2)) algebra is obtained in the limit g — > 0. 

Coalgebra. The defining relations of Q are compatible with the following coalgebra 
structure. The coproduct of the group like elements X G {£7, V, K} is A(X) = X ® X and 
the coproducts of the Chevalley-Serre generators Ej and Fj (J = 1,3) take the standard 
forms. However the coproducts of the fermionic generators E2 and F2 involve an additional 
braiding factor U, which is one of the central charges of the algebra alluded to in the 
previous paragraph, 

A(Ej) = E j ®l + Kr 1 U +s ^®E j , A(Fj) = Fj (g) Kj + U~ Sja <g) Fj. (2.6) 

The coalgebra can be extended to a Hopf algebra. We will give the relevant definitions of 
the antipode and counit later on. 

2.2 Affine Extension 

The infinite dimensional quantum affine algebra Q is the affine extension of Q introduced 
in [35]. The affine extension is obtained by adding an additional node into the Dynkin 
diagram as depicted in Figure 2. The remarkable property of this diagram is that the 
additional fermionic node is a copy of the second node. Therefore, we introduce the affine 
Chevalley-Serre generators {E4, F4, K±} as copies of {E2, F2, K2} and assume that they 
satisfy the same commutation relations as are given in (2.1), (2.3) and (2.4) and also have 
the same coalgebra structure (2.6). Thus, we introduce an additional set of the parameters 
g, a and central charges U, V. We distinguish these two sets by adhering subscripts to them 
arising from the generators to which they are associated, 

g^g k , a^a k , U U k , V -> V k , with k = 2,4. (2.7) 

Next, we need to determine the commutation relations {-E25-F4} and {^4,^2} in such way 
that they would be compatible with the coalgebra structure, 

A({E2,F i }) = {A(E 2 ),A(F A )} and A({E A , F 2 }) = {A(E 4 ), A(F 2 )} . (2.8) 
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Algebra. As a result, we obtain the quantum affine algebra Q [35]. The mixed commu- 
tation relations of it are given by = 1,3) 



{E 2 ,Fi} = -~ga~ l (K A - U 2 U^K^) 



-DAi 



[E j ,F j } = D i 



q - q^ 1 



q ' ' /-A',. 
{E A , F 2 } = ga(K 2 - U 4 U 2 l K^), 

[Ei,Fj} = 0, for i^j, i + j^6. (2.9) 



with the two new constants g and a and the associated supersymmetric Cartan matrix A 
and normalization matrix D given by 



DA 



(+2-1 0-l\ 
-1 0+1 
+1 -2 +1 
\-l 0+1 0/ 



D = diag(l, -1,-1,-1). 



(2.10) 



These are supplemented by the following Serre relations (j = 1,3 and k = 2,4) 

[Ei,Es] = E 2 E 2 = E4E4 = {E 2 , E4} = , 
[F!,F 3 ] = F 2 F 2 = F4F4 = {F 2 , F 4 } = , 
[Ej, [Ej,E k ]] - (q - 2 + q-^EjEkEj = , 
[F 3 , [Fj,Fk]] - {q - 2 + q-^F^F, = . 



(2.11) 



The central charges are related to the quartic Serre relations as 



2,4) 



g k a k {\ - U 2 k V k 2 ) = { [E k , E{\, [E k ,E 3 ]} - (q - 2 + q^E^E^ 
9kcq X {V- 2 - J/," 2 ) = {[F k , Fi], [F k ,F 3 ]} - (q- 2 + q-^F^F^ . 

and the central charges are related with Cartan charges through (k = 2,4) 



v: 



KiK 2 k K 3 . 



(2.12) 



(2.13) 



Coalgebra. The group-like elements X £ {1, ifj, Vfc} (j = 1, 2, 3, 4 and = 2, 4) have 
the coproduct A, the antipode S and the counit e defined in the usual way, 



A(X) = X<g>X, S(X) = X~ 1 , e(X) = l, 



(2.14) 



while the coproducts of the Chevalley-Serre generators are deformed by the central elements 
U k as follows (j = 1,2,3,4), 



A(E j )=E j ®l + K- 1 U2 



Aft) 



l TT + S j,2 TT + S J, 



E 



S(Ej) 
S(F 3 ) 



-JJ~ &j - 2 JJ~ &jA K-F- 



e{Ej) = 0, 



-U. 



+Sj,2 TT +Sj 



Uj u ^FjKj\ e(F j ) = 0. 



(2.15) 



It is important to note that the above coproducts are compatible with all the defining rela- 
tions, including the commutators {E 2 ,F{\ and {E^,F 2 } in (2.9). The opposite coproduct 
is defined as A op = PAP with V being the graded permutation operator. 
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Parameter constraints. In general, the quantum affine algebra Q has seven parameters 
gk,ak,a, g,q (k = 2, 4). A suitable choice of them which lead to an interesting fundamental 
representation was performed in [35]: 

g 2 

92= 94 = 9, a 2 = a 4 a 2 = a, g 2 = - n7 U2 ■ ( 2 - 16 ) 

l-g 2 {q-q L ) 2 

This choice of parameters is also compatible with the bound state representations. Thus 
in this paper we only consider the quantum affine algebra Q, parametrized by four inde- 
pendent parameters g, a, a, q given in the relations above. 

3 Quantum oscillators and representations 

In this section we will provide all the necessary background for constructing the bound 
state S-matrix for the (/-deformed Hubbard model. We will build the bound state repre- 
sentation by introducing g-oscillator formalism linking it to the aforementioned quantum 
affine algebra. 

3.1 q- Oscillators 

We first introduce the notion of (/-oscillators and discuss how to obtain the representations 
of the quantum deformed algebras using (/-oscillators. A concise overview of the (/-oscillators 
and their relation to such representations may be found in [47, 48]. 

Definitions. The (/-oscillator ((/-Heisenberg-Weyl algebra) U q (i) 4 ) is the associative unital 
algebra consisting of the generators {at, a, w, u> -1 } that satisfy the following relations, 

w a ' = q a'w, pa = ato, (3-1) 

ww" 1 = w~ 1 w = 1, a a^ — qa^a = w^ 1 . 

From the defining relations one can see that the element iv^ 1 {a^a — w q Zq-i ) i s central. As 
such, we will set it to zero in the remainder. Then one easily obtains 

ata = ^^, aat = 9 "- g "Y 1 . (3.2) 

We will also need to consider the fermionic version of the (/-oscillator. The above notion is 
extended to include fermionic operators by adjusting the defining relations in the following 
way (we keep the same notation for bosonic and fermionic a, a' for now) 

wa Jf = qa^w, qwa = aw, (3.3) 

ww -1 = w" 1 w = 1, a a^ + qa^a = w. 

In this case, the central element is w 

( a t a _ a gain we set this element to zero, 

resulting in the following identities 

ata = ^^, aJ = qW ~ l -\ lw . (3.4) 

q-q~ l q-q 1 

Of course in the fermionic case the operators a,a^ square to zero. Equation (3.4) implies 
that this only is consistent if w 2 = 1, q 2 . Below we will identify w = q N , where N = 0, 1 is 
the number of fermions making it indeed compatible. 
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Fock space. The g-oscillator algebra can be used to define representations of U q (sl(2)) 
in a very simple way. Let us first build the Fock representation of lA q (§±). For this purpose 
consider a vacuum state |0) such that 

a|0)=0, (3.5) 

then the Fock vector space J- generated by the states of the form 

|n) = (at)»|0) , (3.6) 

is an irreducible module oiU q {\)i). Let us first consider the bosonic ^-oscillators. With 
the help of the defining relations (3.1) and (3.2) one finds that the action of the oscillator 
algebra generators on this module is 

a'\n) = \n + 1), a\n) = [n] q \n — 1), w\n)=q n \n). (3.7) 

This makes it natural to identify w = q N , where N is understood as a number operator. 
Analogously, fermionic generators are found to act as 

a+|n) = \n + 1), a|n) = [2 - n] q \n - 1), w\n) = q n \n). (3.8) 

However, due to the fermionic nature, n can only take the values and 1 and thus the 
identity [2 — n] q = [n] q holds. 

Next consider two copies of bosonic (/-oscillators ai,a\,Wi = q Ni which mutually com- 
mute. Then the Fock space is naturally spanned by vectors of the form 

\m,n) = (a\r(alr\0). (3.9) 

It is easy to see that under the identification 

E = a 2 a u F = a\a 2 , H = N 2 -N U (3.10) 

the Fock space forms an infinite dimensional ZY,j(s[(2))-representation. Moreover, the sub- 
space Tm = span{ | m, M — m) \ m = 0, . . . , M } is an irreducible ZY g (s[(2))-representation 
of dimension M + 1. This can be straightforwardly generalized to sl(ra) and more generally, 
by including fermionic oscillators, this space is extended to the representations of sl(n\m) 
[39]. 

Representations of centrally extended U q (su(2\2)). We will now construct the bound 
state representation for centrally extended U q (su(2\2)) in the g-oscillator language. We need 
to consider two copies of s((2), a bosonic and a fermionic one. Thus we need four sets of q- 
oscillators aj, aj, Wi = q Ni , where the index i = 1, 2 denotes bosonic oscillators and i = 3, 4 
- fermionic ones. Using these we write 

E 1 = alat, F 1 = a\a 2 , H 1 = N 2 -N 1 , (3.11) 

+ + + + iVi + N-i — No — Na 
E 2 = a a\a 2 + b a[a 3 F 2 = c aj,ai + d a3,a 4 , H 2 = -C + — ^— = i , (3.12) 

E 3 = a\a 4 , F 3 = a\a 3 , H 3 = N 4 - N 3 , (3.13) 
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where C is central. It is then straightforward to check that this set of generators forms a 
representation of U q (su(2\2)) on the Fock space when restricting to the subspace of total 
particle number M upon setting 

«« = — rm > oc = — t—, , 0,0= r— p. , ca = , - -= . (3.14) 

[M], [M]g [M]g [M]g V 7 

In the above .8,^3 correspond to the right hand side of the Serre relations (2.12) following 
[10]. As a consequence, the central charges satisfy the shortening condition 



"Ml 2 

Here the ^-numbers are defined as 



[C\*-W=[%]* q - (3-15) 



This way of constructing representations of the centrally extended algebra reminds us of 
the procedure used in, e.g. [30], where long representations were be obtained by twisting 
sl(n\m) in a similar way. 

In the q — > 1 limit the g-oscillators get reduced to regular oscillators and their repre- 
sentations coincide with the superspace formalism introduced in [21]. The identification is 
as follows 

d d 
ai, 2 <-> ■q^, A,2 ">i,2, a 3i4 <-> QQ—, aj 4 <-> 3>4 . (3.17) 

Parameterization and central elements. Introducing V = q c and U as in [35], we 
rewrite (3.14) as 

M _M_ _-, _M_ M_ i 

qM _ g-M ' c = qM _ q -M ' 

which altogether leads to a constraint for [/ and V, 

a 2 (v -2 rj-2 )(1 ^. (v-i^W-g^v- 1 ) (3ig) 

_u ){ i-uv)- {qM _ g _ M)2 . (3.19) 

This constraint agrees with the one in [35] by identifying q — > q M , g — > g/[M] q . The 
explicit parametrization of the labels a, b, c, d shall be given a bit further. 

3.2 AfRne extension 

Next we want to consider the affine extension introduced in [35]. Here we will show that 
our representation allows an affine extension. Analogously to [35] we make the ansatz that 
the affine charges act as copies of E2, F2, H%. I n other words, we set 

E4 = a 4 a\a 2 + 64 a{a 3 , F 4 = c 4 a^ai + d 4 a|,a 4 , H4 = -C 4 + — + ^ 3 N<2 — — . 



(3.20) 
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Checking all of the commutation relations is straightforward. Also, due to the defining 
relations (3.18), the equivalent expressions for the affine representation parameters are 
obtained 

JU _M i _M_ M_ _-, 

_ q 2 V 4 - q 2 y 4 _ q 2 V A - q 2 V 4 

a4d4 ~ q M _ q -M ' ° 4C4 ~ q M _ q -M ' 

a 4 6 4 = ^i(l - U 2 V 2 ), c 4 d 4 = 9 £^(V 4 ~ 2 - C/ 4 " 2 )- (3-21) 

However the commutators between the generators E 2 and E4 and also between F 2 and F 4 
induce relations between a 2 , 04, etc. These are found to be 

a 2 d 4 = W T ^2 - q^Vi 1 ), b 2 c 4 = ^-( q -^U 2 U^V 2 - q^V^ 1 ), 

act , m 1 a/ -, , qd . m 1 m -1 , 

C2&4 = [Mf^^ 2 ' q ~ U * UiV ^> d 2 a A = ^(q-^V^ 1 - q^U^UM), (3.22) 

and agree with [35] upon sending q — > q , 5 — >■ pf] - ; as i n the non-affine case. The tilded 
g, a are not independent but constrained parameters; thus there are 12 constraints for 12 
parameters {a k ,b k ,c k ,d k ,U k ,V k }. 

Hopf algebra and variables. The Hopf algebra structure is just as previously dis- 
cussed in Section 2. Here we will introduce Zhukowksy variables that will parameterize 
the representation labels {a k ,b k ,c k ,d k } and central elements U k ,V k for the bound-state 
representation. Following [35] we choose 

g 2 

92 = 94 = 9, oc 2 = a 4 d 2 = a, g 2 = T t-^. (3.23) 

Note that the powers of q in the expressions above are 1 and not M because g 2 (q — q~ 1 ) 2 
is invariant under the bound state map (g, q) h-> (g/[M] q , q M ), thus these equations are 
identical to the ones for the fundamental representation. 

Also, there is a relation between the central elements of the algebra, 

U 4 = ±£/ 2 _1 , V4 = ±V 2 ~ 1 , (3.24) 

that are called the two-parameter family of the representation [35]. We shall be using the 
plus relation in our calculations. 

The mass-shell constraint (multiplet shortening condition) obtained from the expres- 
sions (3.18) and (3.21) reads as 

(a k d k - q M b k c k )(a k d k - q~ M b k c k ) = 1, (3.25) 

and holds independently for k = 2,4. In terms of the conventional parametrization it 
becomes 

+ + (3,6) 
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where £ = —ig(q — q 1 ). One can further introduce a function £(x) 

CW— X + 1 f_ + J { +1/t . (3.27) 

in terms of which (3.26) becomes q~ M C,(x + ) = q M ((x~). This parametrization leads to 
the following expressions of the labels a k ,b k ,c k ,d k of a 'canonical form': 



n , 9 h _ 9 a k x k -x k 



+ 



M 



/ g 7fc ^ 2 ff , / 9 Vkgq* x k -x k 

^[M] q V k a kg {x + k +£,y VM. ^9 Ik £x+ + 1 ' ^ 



where the central charges are 

tt2 _ 1 X t + £ _ „M x fc £ x fc " t^2 _ 1 " 1 _ „M X k x k + £ 



and the relations between x 2 , 72 and x 4 , 74 are constrained by (3.22) to be 

_i_ _i_ 4. 1 iaj 

x 2 = x , x 4 = — , 72 = 7, 74 = — . (3.30) 

The relation between normalization coefficients 0:2 and 0:4 was given in (3.23). Finally, the 
convenient multiplicative evaluation parameter z for the bound state representation is 

z = q- M C(x + ) = q M COO. (3.31) 

3.3 Summary 

For the convenience of the reader we want to summarize all expressions that will be used 
in the consequent calculations of the bound state S-matrix. We will slightly change the 
notation for parameters related to the fermionic nodes. We rename the representation 
parameters and the central elements of the algebra as 

(o 2 , b 2 , c 2 ,d 2 , U 2 , V 2 ) -> (a, b, c, d, U, V), 

((Z4, 64, C4, c?4, U4, V4) — > (a, b, c, d, U, V), (3.32) 

in order to reserve the subscript position for denoting the momentum dependence, i.e. 
ai := a(pi). in order to reserve the subscript position for discriminating states living in 
different tensor spaces. We will also give some relations that we found to be very useful. 

Explicit representation. The bound state representation is defined as 

\m,n,k,l) = (4r( a lna\) k (4) l \0). (3.33) 

The total number of excitations is k + I + m + n = M. The triple corresponding to the 
bosonic s[(2) is given by 

Hi\m, n, k, I) = (I — k)\m, n, k, I), 

Ei\m, n, k, I) = [k] q \m, n, k — 1, 1 + 1), F\\m, n, k, I) = [l] q \m, n, k + 1, 1 — 1). (3.34) 
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The fermionic part is 

H?\m, n, k, I) = (n — m)\m, n, k, I), 

E^lm, n, k, 1} = \m + 1, n — 1, k, I), F^\m,n, k, I) = \m — l,n + 1, k, I). (3.35) 
The action of the supercharges is given by 

H2\m,n,k,l) = — Ic 1 \m, n, k, I), 

E2\m, n, k,l) = a (— l) m [/] g \m, n + 1, k, I — 1) + b \m — 1, n, k + 1, 1), 
F 2 \m,n,k,l) = c [k] q \m + 1, n, k - 1, 1) + d (-l) m \m,n — l,k,l + 1). (3.36) 

The parameters a,b,c,d are related to the central charges via (3.14). The affine charges 
are defined exactly in the same way, 

H^m, n, k, I) = — Ic — l|m,n,fe,Z), 

E^\m, n, k,l) = a (—l) m [l] q \m, n + 1, k, I — 1) + b \m — 1, n, k + 1, l), 

F 4 \m,n,k,l) = c [k] q \m + 1, n, k - 1, 1) + d (-l) m \m, n - 1, k, I + 1). (3.37) 

The representation labels a, b, c, d are given by 



9 h= I g ax 

[M] q % V[M] 9 7 x- 

M M 

9 7 i Q 2 9 , / g gq^V x 



[M] q aV g{x+ + V [M] q i 91 + 1 



(3.38) 



and the affine parameters d,b,c,d are acquired by replacing V — > V = V x , 7 — > 

a a a 2 and a; 1 * 1 — )• -j-; the corresponding central elements are given by V = (/ , V" = <jr . 



Useful relations. The evaluation parameter z may be expressed explicitly in terms of 
x parametrization as 

• (« - «-')(« - r 1 ) = - [4(- + - x- + ^ - i). (3.39) 
Then using the identity 

Z-t~ 1 = - ( — ^Tu- ( 3 - 4 °) 
)g 

one can further show that it is related to the representation labels (3.38) and their affine 
partners in a very nice way, 

z = -^—(ab-bd), - = ^^(cd-dc), (3.41) 
g a a z g 

while the consistency conditions (3.22) give 

1 _ U2 V 2 1 _ fj2y2 



V 2 -U 2 " v 2 -u 2 



(3.42) 



-13- 



Rational limit. The rational limit is usually obtained by substituting q = 1 + h and 
then finding the h — > limit. Thus by defining the evaluation parameter (3.31) as z = q~ 2u 
we can expand it in series of h as [35] 

z = l-2hu + 0(h 2 ), where u = ^(x + + x~)(l + l/x + x~). (3.43) 

It is noted that the x^ parameters in (3.43) satisfies the leading order of the following 
relation which is stemming from the mass-shell constraint (3.26) in the h — > limit, 

x+ + -!--£-- — = — + 2hMu + 0(h 2 ). (3.44) 
x + x~ g 

In fact, this is consistent with the rational constraint for x^ parameters [27]. Finally, it 
would be important to see how the representation parameters reduce in the rational limit. 
The representation labels (3.38) in the q — > 1 limit reduce to the usual (undeformed) labels 
(a,b,c,d) of [27]. On the other hand, the affine parameters are related to the non-affine 
ones (a, b, c, d) through [35] 



r,i rr ° with m= ° t ™ . t= a -" 

10 I -or x or x \-c d . 



where z is the evaluation parameter given in (3.31), (3.42) and w is defined by 

gv gu 2 -i gg 1/2 u 2 - v 2 

W gq^ 2 V 2 U 2 - 1 gV U 2 - q ' 1 > 

Since the central elements specialize to (U, V) — > {\J~^p, 1) in the limit q — > 1, it is easy to 
see that the matrix relation (3.45) reduces the following simple form, 

MT = T. (3.47) 



4 The S-matrix 



We shall consider the bound state S-matrix which is an intertwining matrix of the tensor 
space furnished by the vectors 

\mi,ni, ki,h) <g) |m 2 , n 2 , k 2 , h)- (4.1) 

Here < mi,ni,wi2) n 2 < 1 an d ki,h,k2,h > denote the numbers of fermionic and 
bosonic excitations respectively with the bound state number Mi being the total number 
of excitations Mj = + m + ki + 1%. Thus the S-matrix is the automorphism of the 
quantum deformed tensor space and is required to be invariant under the coproducts of 
the affine algebra Q, 

SA(J) = A op (J)S, for any J G Q. (4.2) 
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We normalize the S-matrix in such a way that the state 1 0, 0, 0, Mi)® |0, 0, 0, M 2 ) is invariant 
under the scattering. Therefore we will denote the state 

|0) = |0,0,0,Mi)® |0,0,0,M 2 ), (4.3) 

as the vacuum state. 

The invariance under bosonic symmetries AH\ and AH3 requires the total number of 
fermions and the total number of fermions of one type 1 

Nf = mi + m 2 + ni + n 2 + 1l\ + 2l 2 , 

N f3 = m 1 + m 2 + h + h ■ (4.4) 

to be conserved. This conservation divides the space (4.1) into five types of invariant 
subspaces of the S-matrix: 

I \0,l,ki,h)® \0,l,k 2 ,l 2 ), 
lb |l,0,fci,Zi)® |l,0,fc a ,J 2 ), 

II {|0,0,fci,Zi>® |0,l,fc2,i 2 ),|l,l,fci,Ji) ® |0,l,fc 2 ,/ 2 ), 

|0,l,fci,li) ® |0, 0, fc 2 , i 2 >, |0, 1, fci, Zi) ® H,l,fe2,Z 2 )}, 

lib {|0,0,fci,ii> ® \l,0,k 2 ,l 2 ), ® \l,0,k 2 ,l 2 ), 

|l,0,A;i,ii>® |0,0,A; 2 ,Z 2 ),|l,0,A:i,Zi)® \l,l,k 2 ,l 2 )}, 

III {|0, 0, fci, Zi> |0,0,fc 2 ,Z 2 ), 10,0,^1) ® |1, 1, fca, «z>, |1, 1, fci, Zi> ® \0,0,k 2 ,l 2 ), 

|1, l,ki,h) ® |1, 1,^2,^), |0, 1, fci,Zi) ® 1 1,0, A; 2 ,Z 2 ), |1, 0, fci , Zi> ® |0, l,k 2 ,l 2 )}. 

Subspaces I, lb and II, lib are isomorphic, hence we need to find the S-matrix for one 
of the isomorphic subspaces only. In the following we will consider the scattering in the 
subspaces I, II and III only. 

The invariant subspaces differ by the numbers Nfj 3 . By considering the action of the 
algebra charges it is easy to see that the different subspaces are related to each other in 
the way shown in figure 3. 

Finally we want to give a remark on our choice of the basis. The g-oscillator basis we 
are considering is orthogonal, but not orthonormal, 

(m',n',k',l'\m,n,k,l) = * j 5 m ,m' 5 n ,n'5k,k'&i,i' , (4.5) 

where [n]\ = [n] q [n— l] q ■ ■ ■ [l] q is the quantum factorial. We shall choose the normalization 
for the bra vectors to be 

(m,n,k,l\ := \m,n,k,l) j . (4.6) 

1 Note that a bosonic excitation may be interpreted as a combined excitation of two fermions of different 
type. 
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Figure 3. The invariant subspaces of the S-matrix and the algebraic relations between them. 



which helps us to normalize the scalar product to unity and avoid the appearance of 

— 1 /2 

unpleasant numerical factors of the form ([£;]![/]!) in the derivations. The price we 
have to pay for this choice of the basis is that, for real q, the inverse of the S-matrix is 
related to the Hermitian conjugate only up to a basis transformation. For complex q this 
property is not valid even for the fundamental representation [10]. 
For further convenience we introduce these shorthands 

M = M 1 + M 2 , 8M = M 1 -M 2 , K = k 1 + k 2 , 5K = k x - k 2 , 

ki = Mi - hi - 1, Ski = h - h = Mi - 2ki - 1, z\ 2 = zi/z 2 , 5u = u\ — u 2 . (4.7) 

4.1 Scattering in subspace I 

The conserved fermionic numbers (4.4) for the subspace I are Nf = 2K+2 and Nj :i = K+2. 
Thus for the fixed K (0 < K < Mi + M 2 — 2) the dimension of the space is K + 1 and the 
states in this space are defined as 

k 2 ) 1 = |0, 1, k u Mx - h - 1) ® |0, 1, k 2 ,M 2 -k 2 -l). (4.8) 

We start by considering the highest weight state (the state with k\ = k 2 = 0). The 
invariance under AH± and AH% requires it to be an eigenstate of the S-matrix, 

S |0, 0) 1 = &\0,0)\ (4.9) 

Let us compute @. First, we construct the highest weight state by acting with the combi- 
nation AE 2 AE4 on the vacuum state (4.3) (we use the notation Oj = a(pi) etc.) 

Mi ~ ~ T 

AE 2 AE A \0) =q—[M 1 } q [M 2 ] q (a 1 a 2 U 1 V 1 -a 2 a 1 U 1 V 1 )\0,0) 1 . (4.10) 
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This construction let us to rewrite (4.9) as 



S |0, o). . ^ s -^-t |0) 



q— [Mi] g [M 2 ] g (aio 2 UiVi - a 2 a l V ± U h 
A°PE 2 A op E 4 § 

q— [Mi] g [M 2 ]g (oia 2 UxVy - a 2 a x VJJx) 

m 2 -^i a 2 a\ U 2 V 2 — CL1CL2 V 2 U 2 



|0) 



aid 2 — a 2 ai ViC/i 



0,0) : , (4.11) 



where we have used the invariance condition (4.2) when going from the first to the second 
line. Comparing (4.11) with (4.9) we find Q> to be 

s = q ^Mx Q 2 ai U2V2 ~ ai«2 V2U2 = q SM/2 U2V2 xf - z 2 
a x a 2 U1V1 - a 2 di Villi U\V\ x{ - x\ 

In the q — > 1 limit this is the inverse of the result found in [27] due to the interchange of 
A and A op with respect to the ones in [27]. 

Next we define the action of the S-matrix on the subspace I to be 

K 

§\k 1 ,k 2 ) 1 = ^2^ lM \n,K -n) 1 . (4.13) 

n=0 

The strategy for finding coefficients St^n 1 ^ 2 will be based on building the generic state 
\ki,k 2 ) 1 by starting from the highest weight state 1 0, 0) 1 . This allows us to relate ^n 1 ^ 2 
with any k\, k 2 and n to the already known coefficient Q). Thus we need to construct 
k\- and A; 2 -raising operators. We start from inspecting the action of the coproduct of the 
bosonic charge F\ giving 

AF 1 \k u k 2 ) 1 = [h] q q &k2 \ki + l,k 2 ) 1 + [k 2 ] q \h,k 2 + l) 1 , (4.14) 

and 

A^Filfci,**) 1 = [fci] ff |A;i + l,k 2 ) 1 + [k2] q q Skl \h,k2 + l) 1 . (4.15) 

These coproducts do not have the desired properties we want, but are very close. However, 
with the help of E 2 , E% and E± we can construct a new charge with a similar action, 

9 -{E 2 ,[E 4 ,E 3 ]}. (4.16) 



gaa 



We call this new charge 'the affine partner' of the raising charge F%. The action of F\ on 
the state of the form |0, 1, k, I) is 

Fi|0,l,M> =z[l\ q \0,l,k + 1,1-1), (4.17) 
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where we have used (3.39) implicitly 2 . Then it is straightforward to see that the new affine 
raising charge acts on generic states in subspace I as 

AA \h, k 2 ) 1 = zi [h] q \h + 1, k 2 ) 1 + z 2 q 5k i [k 2 ] q \h,k 2 + l) 1 . (4.19) 

And the action of A op F 1 is 

A°PF l \h, k 2 ) 1 = zi q 5k2 [h] q \h + l, k 2 ) 1 + z 2 [k 2 ] q \h, k 2 + l) 1 . (4.20) 

By combining AF\ with AF± we obtain composite operators having the action of the 
desired form - raising k\ and k 2 separately: 



1*1, fe + I) 1 ~ TUT * 1 ^ ~ I*'- fa>'- < 4 - 22 > 
Then by induction we find that the generic state \k±, k 2 ) 1 may be constructed as 
\k k >' - n ^° 1(Zl AFl " ^ AA) nt'o(^A - *2 g'WQ 

Finding J^f 1 '^ 2 is then straightforward. We only need to act with the S-matrix on the 
expression above and sandwich with a frra-vector as 

JT n fcl > fc2 = l ( n ,K-n\ § \k 1 ,k 2 ) 1 . (4.24) 

Performing similar steps as we did in (4.11) and employing the relations 

(A PF 1 -z 2 q Sh A°PF 1 )\n 1: n 2 ) 1 

= [n 2 ] q z 2 (1 - q 5k ^ 5ni ) |m, n 2 + l) 1 + [nt] q {z x q Sn * - z 2 q 6k ') \m + l,n 2 )\ (4.25) 
(ziA^-g^A^FOK^) 1 

= [n 1 ] q z 1 (1 - g ^+^) | ni + l,^) 1 + [n 2 ] q (z ig 5 ^ - z 2 q Sk *)\ni,n 2 + 1)\ (4.26) 



2 For the consistency of the algebra we also give a definition of the 'affine lowering charge' Ei: 

Si = ^{F 2 ,[F 4 ,F 3 }}, Bi|0,l ) fc ) i> = — |0,l,*- 1,1 + 1). (4.18) 

9 2 
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we find the coefficients of the S-matrix in the subspace I to be 



nr=i[^-^nf=r[M 2 -i] 



i 



n-ii[Mi - i] g uT=ii M 2 - j\ q nf=i(^2 - q 1 
fci / 



m=0 



mqk2{n— m) — k\m— k\ 


V 




k 2 




m 




n — m 



m—l 



p=0 p=l+m 
n—m 



x IK 1 



k<2—n—l 
_j2(M 2 -K+n- P )^ Yl ( 

p=— m 



z\2 q 



Mi+2p 



(4.27) 



where = — and the g-binomials are defined as 

Z2 



[b] q \[a-b] q V 



(4.28) 



Apart from the prefactor 2$, this expression only depends on the quotient z\ 2 and on simple 
(/-factors. The expression above has exactly the form that one would expect to obtain by 
an educated guess relying on the one given in [27]. 

Quantum 6j-Symbol. The coefficients J^ 1,fc2 of the bound state S-matrix may be 
regarded as the coefficients which arise in the fusion rule of the irreducible representations 
of U q (su(2)), thus it is expected that the expression (4.27) is related to the quantum 6j- 
symbol, which is the q-deformation of 6j'-symbol and was first introduced in [49]. 

In order to see the relation with the quantum 6j-symbol, we first rewrite (4.27) in 
terms of quantum factorials. This can be done by introducing the notation z\ 2 = q~ 25u 
and using the following identity several times, 



A+B 
I 2 



q 



A-B 



(4.29) 



Secondary, we shift the index of summation m to M\ — 2 — m. After some computation, 
we obtain the following form, 



fel 3 fe2 



- n)(k ^ n+Su+ m [M 2 - k 2 - 1]! [Su + f - 1 - K]\ 



x [hy.ihy.l5u + 5 -f]\[5u- 5 -f-k 2 + n + l]\ 



[Mi 

2 



n-l]\ [8u+ f-l]\ 
-n + 1]! 

V [m + 1]! ([m - Mi + 2 + fei]! [m - Mi + 2 + n]\ [k 2 - n + Mi - 2 



m\ 



m>0 



x [m + Su - 4f + 2]! [Su + % - 1 - m]! [Mi - 2 - m}\ [M - K - 3 - m]!) 1 . 

(4.30) 

where the summation index m runs over the non-negative integers such that all arguments 
of the quantum factorials, which do not include 5u, are non-negative. Finally, replacing the 
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six variables (Mi, M2, ki, k 2 , n, <5u) by the appropriate combinations of (ji, fa, fa, 3a, fa, jo) 
as (see also [27]), 



x (f -2-fc2-<yu), 



•?1 - 2 

J3 = 5 (Mi - 2 - fei - n), je - 5 

we have found that the expression (4.27) obtains a quite elegant form 

yh -73-74+2.75 +j 6 Jh-h+j3)Ui+h-.u-js) \i± + 3 2 ~ ^ ! b'l + k ^jgjj 



k = \{^-l-K + n + 5u), 
j 6 = ±(M 2 -l), 



(4.31) 



x [j's - J4 + is]! [j 3 + j 4 - fa]\ [j 2 - k + faV- [-32 + k + faV- 31 32 33 

34 35 3G 



[1 + h + h + k] ! b'i + k + fa] ! 

(4.32) 

J4 J5 J6 

where we have defined the rescaled quantum 6j-symbol by 

{-l) m [m + l]!(b'i245 - m]\ [ji346 - m]\ b' 2 356 - m]! 

x [m - j'123]! [m - j'345]! ["i - J24e]! ["i - ji56]!) ■ (4.33) 



3i 32 33 
k k fa 



m>0 



Here we have used the bookkeeping notations j abc = j a +3b + 3c and j abcd = ja+3b + 3c + 3d- 
The above expression is related with the quantum 6j-symbol introduced in [49] as 



3i 32 33 
k k fa 



y/2j 3 ~ \y/2fa ~ 1 ("I) 



-ji-jz+Vh+M+js 



x A (ji , j 2 , fa ) A ( ji , fa ,fa)A (j 2 , k , fa ) A (j3 , k , k ) 



3i 32 33 
k k fa 



(4.34) 



where the triangle coefficient A(a, b, c) is defined to be 

^M- c- +t -f + :r^ + "T - 

Rational Limit. In order to find the rational limit of the matrix S£ (4.27) we first use 
the expansion (3.43) for the spectral parameter z. This leads to 

1 



grkite = q 



WU\Ml-i] q VS=i\M2-3\ 



nflilMi " i], tfLiiMi ~ j] q Uli(4lH q + Q M / 2 - l [f ~ l] q ) 
ki 



n—mkiin—ra) — k\rn—k\ 
z 12 1 



m=0 



h 

m 



k 2 
n — m 



' M 2 

ou - — -p 



+ q 



Mi/2 



Mi 



1/2 „Mi/2+p 
12 9 



. Mi 
ou - — - p 



+ q 



A/2/2 



M 2 



m-l / 

xn i 7 V 2/2+p 

p=0 V 

&2— n— 1 

- n 

fci n— m 

x ^ 1_P [M! - p], J] (7 M2 "^ +n ~ p [M 2 - JT + n - p] 

p=l+m p=l 



1 ' 



(4.36) 
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where 5u = u\ — U2- Now we are ready to find q — )■ 1 limit. The g-numbers [x] q coalesce to 
x, thus (4.36) becomes 



U-Li(Mi-i)U k MM2-j) Uf=i(Su + f-l) 

ki / / , \ / , \ m—1 / t- j ^ \ k2—n—l 



X 



£((*]( h )n>+f-.)lf>-^- 

fei n—m \ 

x [J (M 1 -p)]J(M2-K + n-p)\. (4.37) 

p=l+m p=l J 

This result coincides exactly with the expression obtained in [27] 3 

Classical Limit. It is also important to find the classical limit g — > oo of (4.27). This 
limit corresponds to the case 'T(h)' in the analysis of the classical algebra [34], where the 
deformation parameter q is expanded as 

g = l + f + OGT 2 ), (4-38) 
25 

and the parameters become 



x^ 



x ± KM (x + h)(l + l/xh) + 2 
2<7 



ih 

where h = , . (4.39) 



The above expressions are compatible with the constraint (3.26) up to a given order. Since 
£ — > h and x^ — > x in the classical limit, it is easy to see that the evaluation parameter z 
reduces to 

(x + h)(l + l/xh) _ C + D 

h-h-i ~ c-d> [AA[)) 

where elements C and D are the classical limits of U = q D and V = q respectively, and 
are given by 

D = Uz + l)q, C=Uz-l)q, where q = -M H ~ H . (4.41) 



3 The normalization of the evaluation parameter is slightly different in here, Uhcrc = — 2u[27]. 
4 The classical evaluation parameter given in [34] is related with ours as zfg 4 i = (z^'oro) -1 an( l the classical 



'[34] 

parameter is £[34] = —ihh^ 1 (xh CTC + h). 
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nr=i(Mi-;)nf=rr(M 2 -j 



With these preliminaries, we find the classical limit of (4.27) to be 
~ (1 + &d 

ki 



, ki+k 2 M 



m=0 



h 



1 



g Z12 - 1 



fci+n— 2m 



2=1 



-12 



n — m 



9 p=0 

h 



Z\2 ~ 1 



9 ^ 

p=—m 



Z12 ~ 1 



+ 2-(fe 2 (n-m)- fcim-fel) JJ (Afi - p) JJ (M 2 - if + n - p) 

y / p=l+m p=l 

(4.42) 

where S> c i is 0(g^ 1 ) term of 3) in (4.27). Since the binomial coefficients force the index 
m to be m < min{/ci,n}, we will discuss the two possible cases separately. They are the 
n/fci case (off-diagonal sector) and the n = k\ case (diagonal sector). 



Off-diagonal sector. In the case when n is different from k\, it is further classified by two 
more cases - if n is bigger or smaller than k±. Firstly, in the n > k\ case, the leading order 
of (4.42) is 0(g~( n ~ kl ^) with m = k\. Therefore the 0(g~ 1 ) term, which contributes to the 
classical r-matrix, is obtained by setting n = k\ + 1. In this situation, the classical limit of 
(4.42) turns out to be of a simple form, 



ayki,k2 



h z\ 



9 z±- z 2 



-k 2 (M 1 - kt 



(4.43) 



Secondary, in the n < k\ case, the leading order is 0(g ( kl n )) with m = n. Therefore, 
the 0(g~ 1 ) contribution is given by n = k\ — 1. In this case the amplitude becomes 



h z 2 



-k 1 {M 2 -k 2 



I) • 



9 zi - z 2 

The other matrix elements do not contribute to the classical r-matrix. 



(4.44) 



Diagonal sector. This is the n = k\ case and it needs a more elaborate treatment in 
comparison with the off-diagonal sector. In this case the leading order in (4.42) is 0(1) 
with m = k\ = n. Thus the classical limit turns out to be 



fcl 



1 + @d - ^{k\ + k\) + -— — 
2g gzi- z 2 



(M , 

i=i v 



fei-i 



k 2 — k\ — 1 



, ^ fz l M 2 -z 2 M 1 \ 

+ E( 2 + ^P)+ E 

p=o v 7 p=-k 



z x M x - z 2 M 2 



+ zip 



(4.45) 
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Full Rational Limit. It is noted that the classical limit still depends on the deformation 
parameter h. This allows us to take h — > limit further, which corresponds to the case 
"R(full)" in the analysis of [34]. In this limit, the classical evaluation parameter (4.40) 
reads, 



2~ 1- -u + 0(h 2 ), with 
9 



u = x + -. (4.46) 

x 



Then the off-diagonal elements of the classical r-matrix (4.43) and (4.44) turns out to be 



On the other hand, the diagonal elements (4.45) reduce to 



-fcl,fc2 



(4.47) 



fci+A:2 fci— 1 k2 — ki — 1 

E (¥-0 + E(-^+p)+ E 

Z=l p=0 P=— fel 



(4.48) 



The above expressions (4.47) and (4.48) agree with the classical limits of rational case [27]. 

4.2 Scattering in subspace II 

The S-matrix in the subspace II is defined to be 

K 4 



S = EEK (^n 1,fc2 ) ■ , 

n=0 j'=l 



(4.49) 



and the standard AN + 2-dimensional basis is 



|fci,*2>? = |0,l,fci,Mi - h - 1) ® |0,0,£; 2 ,M 2 - fc 2 ), 
Ifci.fe)? = |0,0,fci,Mi - fci) ® |0,1,A!2,M2 - Ai2- 1), 
1^1,^2)3 = |0, 1, k\,Mi - h- 1) ® |l,l,fc2 - 1,M 2 -fc 2 - 1), 
*2>4 = |1, 1, fci - 1, Mi - fcx - 1) ® |0, 1, fc 2 , M 2 - k 2 - 1). 



(4.50) 



We shall express the coefficients (^n 1,k2 ){ in terms of already known J^f 1 '^ 2 with the help 
of the charges AZ? 2 and AE4 that relate the states in the subspace II to the states in 
subspace I: 

AE 2 \h, k 2 )f = Qj(ki, k 2 ) \h,k 2 )\ AE A \h, k 2 )f = Qj(k u k 2 ) \kxMf- (4.51) 

The coefficients Qj(ki,k 2 ), Qj(ki,k 2 ) and their partners for A op E 2 and A op £ , 4 are spelled 
out in the Appendix A.l. 
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The strategy of finding ^n 1 ^ 2 is the following. We start by considering the matrix 
element 

4 K 

1 (n,K-n\A op E 2 8\k 1 ,k 2 )f = Y,Y1 J (n,K- n\ A op E 2 \m,K- m)f {&t M )\ 

j = l 771=0 

4 K 

= HY^ l {n,K-n\m,K-m) 1 Q° p (m,K-m) (&t M )\ 

j = l 771=0 
4 

= ^Qf(n,K- n ) {&t M )\- (4-52) 
i=i 

Next, using the invariance of the S-matrix A° V E 2 S = §Ai?2, we rewrite (4.52) as 
J (n, K- n\^AE 2 \k h k 2 )f = : (n,K- n\S\h, k 2 ) 1 Qi(k u k 2 ) 

N 

771=0 

= ^ lM Q l (h,k 2 ). (4.53) 

Likewise we get a similar set of relations by considering the charge E±. These relations can 
be conveniently summarized in terms of a matrix equation 



Qf{n,K-n) Q° p (n,K-n) Q°/(n,K-n) Qf(n,K-n)\ hl 



A- 



Qf(n,iT-n) Q2 P (n,iT-7i) Q° p (n,K-n) Q?(n,K 



n 



( Qi( k i,k 2 ) Q 2 {h,k 2 ) Q 3 {ki,k 2 ) Qi(ki,k 2 )\ 
\Qi(ki,k 2 ) Q2(h,k 2 ) Q 3 (ki,k 2 ) Q4(ki,k 2 )l ' 



(4.54) 



giving a total number of 8 constraints. However, there is a further need of 8 more con- 
straints. These can be obtained by considering a composite operator 



E 2 = e (ex F\F^F 2 + e 2 F 1 F 3 F 2 + e 3 F 3 F 2 F\j , (4.55) 

where 

eo = ^^^(/^i - q 2K+2 z 2 )-\ ei = (q- g" 1 ), 

e2 = q ^( q -*-™ Zl - q *- M z 2 ), e 3 = - q M ^{ q ~^ K z x - q^ M z 2 ), (4.56) 

and its anine partner E4. These operators act on the states in the subspace II as 

AE 2 \h,k 2 )f =Z i (k 1 , k^hM) 1 + Z+(h, k 2 )\h + 1, k 2 - l) 1 

+ Zr(k 1 ,k 2 )\k 1 -l,k 2 + l)\ (4.57) 
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giving 

4 K 

\n,K-n\ A op E 2 S \k 1: k 2 )f = Y,Y1 \ n i K ~ n \ A ° P ^2 \m,K-m)f {^ tM )\ 

j=l m=0 

3=1 V 

+ Z~^(n + l,K-n-l) (^t 2 )l) ■ (4-58) 

The coefficients (4.56) are chosen in a such way that the 'non-diagonal' part of this relation 
is vanishing, Z^~' op (n — 1,K— n + 1) = ZJ' op (n + 1,K— n — 1) = 0. Therefore the only 
surviving part of (4.58) is 

4 

\n,K-n\ A°PE 2 S \k u k 2 )f = ^ Z°f (n, K - n) {®}* M )[. (4.59) 

3=1 

This results in the following matrix equation for Z° p (n, K — n): 

(z op (n,K- n) Z° p {n, K- n) Z° p (n,K- n) Z° p (n, K- n)) 
= (z^hM) Z 2 (h,k 2 ) Z 3 (h,k 2 ) Z A (h,k 2 )) %t M 

+ (z+(h,k 2 ) o z-(h,k 2 ) o) %t+ XM - x 

+ (0 Z£(kiM) Z^(h,k 2 )) (4-60) 

plus a similar set of equations arising from the affine charge E\. Both sets can further be 
united into a compact matrix form 

4#j*i.*a = B%t M + B + ^ 1+1M " 1 + B-3£*-- 1 ' ki+x , (4.61) 

which multiplied from the left by A~ l defines all coefficients of ^f 1,fe2 in terms of already 
known %t' k \ %~n 1±1MTl . The explicit expressions of matrices A, A 1 , B, B^, their 
q —> 1 limit and the coefficients Zi(ki,k 2 ), Z° p (n, K — n) and their affine partners are 
spelled out the Appendix A.l. 

To finalize we want to note that not all of the constraints in (4.60) are linearly inde- 
pendent. The set of independent constraints is chosen in such way that the inverse matrix 
A~ l would exist. 

4.3 Scattering in subspace III 

We will compute the S-matrix components in the subspace III in a very similar way as 
we did in the previous section for the scattering in subspace II. We start by defining the 
S-matrix for the subspace III as 

K 6 

§ \h,k 2 )f l =^2^2\n,K- n)f {^t M t , (4-62) 

n=0 j=l 
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where the standard basis for the 6iV-dimensional vector space is 





fa 




\o,o,fa 


Mi 


-fei) 


® 


0,0,fe 2 ,M 2 - fa), 


h 


fa 




\o,o,fa 


Mi 


-fei) 


® 


l,l,fa-l,M 2 -fa-l), 


fa 


fa 


>§" = 


\i,i,fa 


-1, 


Mi - 


fei 


- 1) ® |0,0,fe 2 ,M 2 - fa), 


fa 


fa 


>T = 


|i,i,fei 


-1, 


Mi - 


fei 


- 1) ® |l,l,fe 2 - - fa 


fa 


fa 




|l,0,fci 


-1, 


Mi - 


fei 


) ® |0,l,fe 2 ,M 2 -fa- 1), 


fei 


fa 


> l i l = 


|0,l,fci 


Mi 




- 1 


) ® |l,0,fe 2 -1,M 2 - fa). 



(4.63) 

Next we shall employ the same strategy as before. We perform the same steps as in 
(4.52) and (4.53) only with A op E 2 , giving 



f(n,K- n\ A°*>E 2 S \fa, fa)f = ^(G op (n, K - n))]^' 1 *) 1 ., 



l=i 
i 



f{n,K-n\SAE 2 \fa,k 2 )f I = 



m=l 



yjG(k U k 2 ))™, 



(4.64) 



where are the matrix representations of the charges A^ op ^E 2 . Once again these equa- 
tions (together with the affine ones coming from E4) do not provide enough constraints to 
define the matrix ifjf 1 '^ 2 uniquely, and we need additional constraints. They are obtained 
with the help of A ( - op )(F 3 F 2 ), namely 

6 

f(n - 9 U K-n + 0i-l\ A op (F 3 F 2 ) § \fa, fa) 1 / 1 = ^(# op (n, n - K))\{^ M )\ , 

1=1 
4 

fin - U K-n + 6 t - 1| § A(F 3 F 2 ) \fa, fa) 111 = £ (W^JjHifa, fa))J , (4.65) 



m=l 



where 6{ is defined by 9i = (1 — (— 1)*)/2 and H^ op ^ is the matrix representation of 
A^ op ^ (F 3 F 2 ). Here we have also introduced & 2 as 



These equations may be written in a compact way using matrix notation 

G op (n,K- n) ^ M = G{fa, fa), 

H op {n,K- n) = H(fa,fa). 



(4.66) 



(4.67) 



The explicit realization of the matrices in the expressions above are spelled out in the 
Appendix A. 2. 

Similarly as in the previous case, not all rows and columns of G^) and arc 
linearly independent, thus we have to select the independent ones only. Therefore by 
taking the following linear combinations, 



G 



(op) 



^K—n — Tf- 



¥ (a 2 G (op) - a 2 G {op) ) and H {op) = c 2 VxH^ - c 2 V^H 



K°p) 



(4.68) 
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where the tilded matrices are the affine counterparts and selecting the first three rows of 
each, we are able to combine them into the non-singular quadratic matrix A (6 x 6) and 
the rectangular matrix B (8 x 6) as follows (j = 1, • • • ,6), 



!=1 ' 2 ' 3 ' and (*)-=( ( 5' 4 !=1 ' W - (4.69) 
. {H )f > ' = 4,5,6, " \(ff)J- 4 , i = 5,6,7,8. 

This approach let us to rewrite the constraints (4.67) in terms of a single matrix relation 

A = &* lM B giving = A~ l &* lM B . (4.70) 

This relation let us to obtain any matrix element (^f 1,fc2 )*- of the scattering in the subspace 
III. Here we have also introduced the block diagonal matrix ty*' 1 (6 x 8) as 

'(^ n M )J, i = 1,2,3, and j = 1,2,3,4, 
{&n l )iZ 3 4 , * = 4,5,6, and j = 5,6,7,8, (4.71) 



, the rest 



The explicit form of matrices A, A 1 , B and their q — > 1 limit are given in Appendix A. 2. 
5 Special cases of the S-matrix 

In this section we consider the reduction of the S-matrix in the case when one or both 
factors of the tensor space (4.1) are transforming in the fundamental representation. 

5.1 Fundamental S-matrix 

As a most simple case of the derivations presented in section 4, we want to compute the 
fundamental S-matrix found in [10]. The fundamental representation is defined by setting 
Mi = M2 = 1 and the corresponding S-matrix is 16 x 16 - dimensional. In order to make 
the comparison with [10] more explicit, let us denote 

a{ >2 = (j) 1 ' 2 , and a\ A = ip 1 ' 2 . (5.1) 

Then, starting with the subspaces I and lb, we find 

SIV^V") = @\i> a 4> a ), (5.2) 

where *3l is given by (4.12). Further, due to our normalization 

§\cj) a cp a ) = \4> a 4> a ). (5.3) 

Here we would like to remark that our normalization differs from [10] where the S-matrix 
is normalized such that §1^"^") = —\ip a ip a )- In other words, the quantities given here 
need to be divided by an additional factor of '3. 

Next we proceed to the subspaces II and lib. For the subspace II (and analogously 
for lib) the parameters fei, /c2, n indexing the matrix W can take the values and 1, but 
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fortunately, we find that <3f is the same for both of these values. Next it is easy to observe 
that the matrices A (A. 4) and B (A. 5) get reduced to the upper left 2x2 blocks 

while the matrices B + and B~ do not contribute at all. This gives the following solution 
of (4.61) 

( v ^(a 2 ai UfVf-ai&UiVi) aiai(l-i7|V 2 2 ) 
Ui Vi (0,20,1— ai a.2 U2V2) a-iazU'zVg— «2ai 

a 2 a2U2(UfV?-l)y2 {a2a 1 -a 1 a2)U2V2 
UiVi(o2ai— ai^t/f V 2 2 ) v / g(a2fii -ai^U'^V^) 



7i U2V2 £2 

72 x i ~ x i 



1 1 / 2 U 2 V 2 ^- = 

1 2 2 X+-X~ 72 C^lVl x +_ 



(5.5) 



c~ q 1 / 2 U 1 V 1 x+- 



Then the corresponding explicit form of the fundamental S-matrix acting on the inequiva- 
lent states is 

%V> = q l ' 2 u 2 v 2 %^\r<?) + ^ , 

72 tvi Vi — x 1 q i ' z UiVi X2 - x l 

Finally we turn to the subspace III which is four dimensional in this case. Analogously 
to our strategy presented section 4.2, we inspect the action of AE2 and AE4 obtaining 

Ai^l.O)? 1 = ^a 2 |l,0)f, A£? 2 |l,0>| n = 6x11,0)° 

A£; 2 |0,l)i n = a 1 |0,0) I 1 I , AE 2 \0,l) 1 i I = -U 1 V 1 ^qb 2 \0,0)f, (5.7) 

plus similar expressions for E4. For completeness, let us spell out the opposite coproduct 
as well 

A^I^O)? 1 = a 2 \l,0)f, A°^ 2 |1,0)P = b 1 U 2 V 2 Vq\l,0) l l 

A o ^ 2 |0,l)i n = oi^|0, 0)?, A^ICl)™ = -6210,0)?. (5.8) 

The equation (4.70) in this case becomes 



w ,u )s , (5.9) 



\a 2 b ly /qU 2 V 2 J (&r)l J \ 

the explicit solution of which is 

(/ av>lfi\l r avl,0\5\ ( 0-~ x 2 x i)( x t ~ x t) x i a ( x i ~ x t)( x 2 ~ x t )( x i~ x t) 

)l )\ J _ [ (l-x^x^)(x7-x%) qxf ^g^Vm^i^ x 2 - 1 )( x i- x t) 

(9fh°)\ \ 7i72(4-a:+) x~ (l-x~ x+)(x+ -x+) U 2 y 2 x 2 

' \ U 2 V 2 a(l-x-x-)(x+-x-) q 3 / 2 x+ (l-x~ x~ ){x~ -x+) U1V1 ~qx~^ 

(5.10) 

The remaining matrix elements are then easily deduced from similar derivations. These 
results are in agreement with [10]. For a complete list of all the scattering elements we 
refer to the Appendix B.l. 
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5.2 The S-matrix S Q1 

In this section we will derive the S-matrix describing the scattering of an arbitrary bound 
state with a fundamental one, Sqi- Once again, we will follow the derivations performed 
in section 4 step by step. First, by setting M 2 = 1, we find that the states in subspaces I 
and lb scatter almost trivially 

SIM) 1 = ®\k,o) 1 . (5-11) 

However the scattering in the subspace II does not get simplified that much. Nevertheless, 
for fixed k\ + k%, the corresponding vector space gets restricted to 

{\h,0)Y, lfci-1,1)! 1 , \h,0)f, \h,0)U}. (5.12) 

This is because the states \k±, ^2)3 have M 2 > 2 and thus they are not present. By reducing 
our general expressions to accommodate these 4 states, we are lead to 16 inequivalent 
scattering elements, however we found 2 of them to be vanishing. The rest may be casted 
in quite compact form as 

SIM)! 1 = (% k '°)\\k,0)Y + (^°)\\k-l,l)Y + (% k '°)l\k,0)l l + (<'°)f|fc,0)? , 

s|*-i, i)! 1 = (^ 1,1 )il^o)! 1 + (sf-^Wk-i, 1)! 1 + (^ 1,1 )?l^o) I 2 I + (^ 1 ' 1 )!|fe,o)i I , 

S|fc, 0>? = (% M )^,0) I 1 I + (W^Hk-l, l)! 1 + (% k '°) 2 2 \k,0) 2 l , 

S|fc s 0>? = (% k '°)l\k,0)Y + (9f'°)i|fc-l, l)! 1 + (% k,0 )i\k,0)% . (5.13) 

The explicit expressions of the coefficients above are given in Appendix B.2. Upon setting 
Mi = 1 the coefficients with indices 1 and 2 reduce to the ones of the fundamental S-matrix 
(5.6) derived previously. 

The scattering in the subspace III simplifies considerably. It is easy to see, that the 
states Ifcijfe)™ n eed not to be considered. Thus we are led to the reduced case of our 
general expressions for subspace III that involve the states (5.12) and 

{|M)! n , \k,0)f l , \k,0)f\ |fc-l,l)i n , \k-l,l)^, \k-l,l)f 1 } (5.14) 

only. However, there is a more straightforward way to obtain the S-matrix in this particular 
case. 

There are 36 scattering coefficients in subspace III that need to be determined, but 
not all of them are independent. Firstly we can relate the half of them to the other half 
by considering the identity 

AE 3 \k - l.O) 1 = M)™ + Q~'\k ~ 1,1)™, (5.15) 

giving 

S |fc - 1, 1) HI =®(\k, 0)f l +q\k-l, -qS\k, 0)f l . (5.16) 
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Subsequently we can express the states \k — 1, l)* 11 , \k — 1, 1)™ as follows 

I^.O)? 1 = |Ar- 1, l>i IT , 



fc,0>J I = |fc-l,l>? 1 . (5.17) 



AF 1 A£ 1 -#-iyM-fc], |; _ MlII _ | ,_ 1 m 



The explicit constraints that follow from these identities are listed in the Appendix B.2. 

Then instead of reducing the general expression of the matrix we follow its deriva- 
tion path. By considering the action of the charges F 2 and F4 on the subspace II states we 
are able to find simple expressions that relate subspaces III to subspace II as 

m n = W-^ |M)?| |M)T n = W^ M)f] 

cid 2 UiV2 - cid 2 UiV2 d\diU\V 2 - d 2 diUiV2 

l K \q c\diV\V 2 - c\diU\Vi 

This approach let us to find the expressions of the matrix elements of 2£ in terms of the 
matrix elements of <3f for this particular case in quite an easy way. The explicit expressions 
are once again given in the Appendix B.2. 



6 Discussion and outlook 

In this work we have constructed the supersymmetric short representations of the quan- 
tum affine algebra Q based on the centrally extended su(2|2) algebra by making use of the 
quantum oscillator algebra. These representations are of great importance as they accom- 
modate the bound states of the model. We found that the bound state representations 
of the affine extension show a lot of similarities with the fundamental one constructed in 
[35]. In particular, we found that the affine central elements are inverse to their non-affine 
partners, exactly as for the fundamental representation. Moreover, the parameterization 
can be derived from the fundamental one by the simple map (q,g) — > (q ' ,g/[M] q ). 

The affine extension plays a key role in the construction of the bound state S-matrix. 
Indeed, the affine generators E4 and F± are crucial in constructing the blocks and <3f . 
In other words, the bound state S-matrix is uniquely fixed up to the overall dressing phase 
by requiring invar iance under the affine algebra Q. 

We have also spelled out the explicit coefficients of the S-matrix when one of the spaces 
corresponds to the fundamental representation. And in particular, we have checked that 
our formalism correctly reproduces the fundamental S-matrix found in [10]. Furthermore, 
our results are in a very good agreement with those of [27], where a similar derivation based 
on the Yangian symmetry related to the same underlying Lie superalgebra was performed. 
More precisely, the S-matrix we have obtained in the q — > 1 limit for the subspace I reduces 
to the one found in [27]. However we can not make a direct comparison for subspaces II 
and III as the intermediate expressions are different. This is due to the fact that affine 
rather than Yangian generators are used. Nevertheless, the expressions we have obtained 
in this work are of more symmetric form than those of [27]. This is an expected result, 
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as the deformed quantum affine algebra itself is of more symmetric form than its Yangian 
limit. 

We have not checked the Yang-Baxter equation in full generality due to this being 
extremely challenging from the technical side. However, we have performed a series of 
checks for a wide variety of states using numerical computations and found that it was 
perfectly satisfied. A more detailed discussion on this point is given in Appendix C. 

In order to complete the investigations concerning the S-matrix it would be interest- 
ing to consider the crossing symmetry and the corresponding solutions for a g-deformed 
dressing phase, which at the moment are not known. 

A particularly interesting direction for future research would be to study the repre- 
sentations and their S-matrices for q being a root of unity. It is well known that the 
representation theory for these values of q differs from the one for real q. Due to the bound 
state map being of the form q — > q M , it is not difficult to see that there appears to be some 
intrinsic periodicity to these representations. One could hope, for example in the context 
of the thermodynamic Bethe ansatz, that this would result in a finite number of bound 
states, leading potentially to some useful insights in this area. 

A different topic related to this, would be to investigate the algebraic Bethe ansatz 
and the bound state transfer matrices. This could perhaps be used to find a g-deformed 
version of the T-system. 

One more possible direction of investigations is to consider the boundary conditions 
and boundary scattering for the deformed Hubbard Chain. A good starting point for 
this approach would be to consider the boundary conditions equivalent to the ones of the 
Y = and Z = giant gravitons in the framework of the AdS/CFT correspondence [50]. 
We expect some sort of deformed (twisted) coideal subalgebra of Q to be governing the 
boundary scattering of the aforementioned type that in the rational limit would reproduce 
the twisted Yangian algebras constructed in [51-53]. 

Other open questions include the search of the complete algebraic R-matrix and a de- 
tailed investigation of its classical limit along the lines of [34]. It would also be interesting 
to extend the classical limit to the next order. For the undeformed case it was found that 
this order coincides with the square of the classical r-matrix [54]. 
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A Elements of the S-matrix 



In this Appendix we have spelled out various coefficients and matrices that have been 
heavily used in the intermediate steps in deriving the final expressions of the S-matrix for 
the subspaces II and III. 

A.l Subspace II 

The coefficients for the charge AE 2 in (4.51) are 

Qi(h, k 2 ) = -^/ 2 - fc i a 2 UxVx [h + l] q , 

Qz{kiM) = -q Ml/2 ~ kl b2U 1 V 1 , 
Similarly, the coefficients for the charge A op E 2 are 

QT(ki, k 2 ) = -a 2 [k 2 + l]q, Q° 2 {ki, k 2 ) 



Q 2 (kt,k 2 ) = ai [ki + l] q , 
Qi{ki,k 2 ) = bi. 



M 2 /2-k 2 



ax U 2 V 2 [ki + 1] 



q 1 



Q 3 (ki,k 2 ) — —b 2 , 



Qf( 



<i,k 2 ) = q hhl2 ~ k2 hU 2 V 2 . 



(A.l) 



(A.2) 



By replacing a,b — >■ a,b and U, V — > U,V, one obtains Qi(ki,k 2 ) and Q° p {ki,k 2 ) related 
to the affine charge E±. 



The coefficients in (4.60) are 



and 



Z?{n,K- 
Z°/(n,K- 

1 

Z 1 (k 1 ,k 2 ) 
Z 2 (ki,k 2 ) 
Z 3 (ki,k 2 ) 
Z 4 (ki,k 2 ) 



n) 
n) 



c 2 Vt [M 2 -K + n] q , 
d 2 Viq~ M \ 



Z° p {n,K- 
Z° p {n,K- 



n 



c 1 U 2 [n-M l ] q q n - 
-d x U 2 f 



Hi ' 



(A.3) 



c 2 Ui [k 2 + l] q ^ Ml/2 _ kl+ M 2 



q M z 12 - q 2 ( K + 1 ) 

Z12 ClV 2 [h +l]g 



,2(i 



Z12 - q 



SM( n 2(n-ki) 



1 



q M z 12 



q -SM/2+2l q 2n Z2i 



d 2 U x 



,M 



-,M 



z X 2 ~ q 2{K ^ 
z\2 d\V 2 



qMl /2- kl ( 2n zi2 



- g 5M (g 2 ( n + fc 2) 

q M/2(n-h) _ 



„.2K- 



-2k 2 +5M 



2k 2 +5M 



1) 



Z12 



7 2(if+l) 



M/2+2 ( 2n 



Z21 - q 



-Ml 2(n+k 2 ) _ 2K\ 



- q 

,2k 2 +5M 
2k 2 +5M 



The matrices in (4.61) are defined as 



A 



(QT(n 


K- 


n) 


Q° 2 (n 


K — 


n) 


Q° 3 p {n 


K- 


n) Q op (n 


K- 


n)\ 


QT(n 


K — 


n) 


Q° 2 {n 


K- 


n) 


Q° 3 p {n 


K- 


n) Q° p (n 


K — 


n) 




K — 


n ) 


Z° p (n, 


K— 


n) 


Zj p (n, 


K — 


n) Zl P (n, 


K — 


n) 


\Z° p (n, 


K- 


n ) 


Z 2 p (n, 


K- 


n) 


Z° p (n, 


K- 


n) Z2 P (n, 


K- 


n)J 



(A.4) 



B 



fQi(h,k 2 ) Q 2 (h,k 2 ) Q3(ki,k 2 ) Q4(ki,k 2 )\ 
Qi(ki,k 2 ) Q 2 (fa,k 2 ) Qz{ki,k 2 ) Qi(h,k 2 
Zi(k 1 ,k 2 ) Z 2 (k 1 ,k 2 ) Z 3 (ki,k 2 ) Z i (k 1 ,k 2 

\Zi(ki,k 2 ) Z 2 (k!,k 2 ) Z 3 (ki,k 2 ) Zi(ki,k 2 )J 



(A.5) 
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and 



B~ 



( o o\ 


Z+(kx,k 2 ) Z+(k u k 2 ) 
\z+(kx,k 2 ) Z+{k x M) 0/ 



/o \ 


Z^{k\,k 2 ) Z^(ki,k 2 
\o Z^(h,k 2 ) Zl{k u k 2 )J 



The latter two have a quite compact explicit form 



B + 



q 



l+k 1 -k 2 - 



Mi 



2- q Kh+2k 2zi2 _ q M 2 +2(n+l) 



(q-q- 1 ) 



l M z X2 - q^ K ^) 



B~ 



„l-ki+^f q M 2 +2n Zi2 _ Afi+2(fe 2 +l) 



(q 



j M z 12 - q 2 ( K +V 



o o\ 

0^ 
-c 2 Ui[k 2 ] q d 2 Ux 
\-c 2 Ui[k 2 ] q d 2 Ux 0/ 

/o \ 

0^ 
-ciFaNg o diy 2 

\o -cxv^kxig o Jiy 2 y 



(A.6) 



(A.7) 



(A.8) 



The inverse of A has a very complex form, however it can be decomposed intro three quite 
compact matrices as A -1 = CVD, where 



C 



( 



z\ 2 b 2 







z\ 2 ab 2 



[M 2 -K+n] q " [M 2 -K+n] q 

q K --r- n ab 1 U 2 V 2 n 





M 2 





-z\ 2 a 2 




[n-MiJg 


M 2 



p7^ 







n 5aiC/ 2 ^2 



ai 



C/2V2 



(A.9) 



D = diag 



M 2 

q 2 



gaa 2 z 2 ViV 2 a V1V2 



(A.10) 



r 



here 



1 

TP 



Z12 



C/ z £ 2 -V 2 + 



v z -u z 



iiu z 



-V, 



\{v z -u z e) v z m z i 



~ Vz+ z 12 

v z -u z 



iUzi ~V Z J 



(A.ll) 



w = v z v z - u z u z f, u z 

plus similar expressions for V z . 



Z12 - Uful 



Z12 - Ulul 



(A.12) 
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Rational limit. The matrices B + (A. 7) and B~ (A. 8) in the q ->■ 1 + h (h 

become 



B~ 



2hk 



5u-^-k 2 + n + l 



5u 



M 



+ K + 1 



B 



- Su + ^ + k 2 - n + 1 
5u - f + K + 1 



/ o o\ 


-k 2 c 2 /Ux da/E/i 
^202^1/00 —62^1/0" J 

/o \ 


— k\c\ di 
\^0 k\a\ja.a. — b\/aa J 



The matrices ^4 (A. 4) and -B (A. 5) in the q — > 1 limit become 



.4 



B 



-(M 2 -K+n)g 2 ~t2 (Mi-ra)gil/ 2 7i 2 _ 2 - - L ^ » 

72^2 Tl^l 
i(M 2 -X+n)qg272 i(A/i-n)ggi7i ia&gg^j-^) iaagifo^-sj') 

a;+ U2~x^~ 72 ^271 

i(M2-K+n)g 2 -y2 iJMj-n)grn mix? -x%) igi(x^-x^) 

ax% aU2X J [ 72 ^271 

(M 2 --R'+n)g272 (Mi-n)g 1 t/ 2 7i £2(a£-£^) gi?7 2 (a:~-3;+) 

Qa a " 072X2 0712:^ 

ag 2 Ui(x^ -x^) agi(a:~-x+) ^ 

72^2" 7i a; r 
i(M2-k2)ag2-f2 t(Mi— fci)agi7i ia«g2(:r2 -^2") iaagi^r -1 !") 

(7ix+ 1+ ^172 71 

i(M 2 -fc 2 )g272 t(Mi-fci)gi7i ^(x^-^) tgi(xj'-gj') 



'(M 2 -k 2 )g 2 U 1 'y2 {Mi — fci)gi7i 



aUix^ 



11 



V 



(M2-k 2 )g2Ui72 (Afi-fci)gi7i g2t/i(x 2 -sj") gi(x 1 



0) limit 



(A43) 



(A44) 



(A.15) 



(A46) 



The notation used in here is gi = y -^f and Ui = 

It might seem that the matrices B + and B~~ do not contribute in the q — > 1 limit 
as they are of order 0(/i), however the combinations A~ 1 B + and A~ l B~ in (4.61) are of 
order 0(1), thus are defined correctly. We do not spell out the explicit expression of A" 1 
in the q — > 1 limit as it is quite sizy and also not much illuminative. 

A. 2 Subspace III 

The coefficients' matrices in the expressions (4.67) 

G op {n, K- n) ^ M = G{h,k 2 ), 
H°P(n,K- n) 2?t M = H(k u k 2 ), 
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arc 



G op 



-K+n 



[Mi—n] q Ol 



M 2 



K+n 



U2V2 
[Mo-K+n} q a,2 






- x+n [Mi^] q ai 



U2V2 




U2V2 















M2_ 



■K+n, 



~b 2 





M 2 



[M2-A"+n] 9 a2 



C/2V2 
&2 



C/2V2 

[M 2 -K+n] q a 2 



^--K+n 



[Mi—n] q ai 



G 



[Afi— fci] g oi 



Mi , Mi , 



62 



C/iVi 






C/1V1 
[Ml-fei] g ai 











61 



-C/2V2 





61 








Mi 



/ 

(A.17) 
62 \ 



-C/1V1 




Mi 



J l 1 

C/iVi C/iVi 



fc l[M 2 — fc2]a2 

C/1V1 



and 



■\q£l_ 



U 2 



'[K-n] q c 2 q n 





Mi 



C/ 2 



-<h 



Vi 




-Vi 

C/2 











Mi 



-[Afi— fci],oi 



Mi 







(A.18) 



Vi 





" [K-n] q c 2 



( a k ^ 



H 



" [fcl]gCl 



Vi 



Mo 



-Vi 



Vi 




f/2 



Wg£l 
C/2 



V2 

[fc 2 ] 9 c 2 c7i 







-d 2 c7i 

Mo 





-Va 




[^2]qC 2 (7l 





Mo 



"di 



-V 2 





[A; 2 ] g c 2 t7i 




(A.19) 



7 *2" 



V 2 



fc - A/2 
2 ~2~ [fcl]gCl 



v 2 



(A.20) 



Their affine counterparts G, G op and i?, ff op are obtained by the replacing non-afhne (or 
affine) parameter to affine (or non-affme) ones. The matrix <3f n ' is a slightly modified 
version of ' ki , 



W kl ' k2 





l,fca 


){ 




- X )\ 




l,k 2 


)l 








l,fea 


)? 








IM 


)§ 








l,fea 


)? 


^<yyki,k 2 - 






IM 


)I 


(ayk ly k 2 - 






l,fea 


)f 






{Wn 1 - 


IM 


)i 


[& kx,k ' 2 ' 





(A.21) 
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The coefficient matrices in (4.70), A 2f r 



ki,k 2 



&n' l B, are 



.4 



/ [M!-n]A 3 
1 U 2 V 2 




[n} q A'2 



V 



U 2 Vl 








-Q2Z2 
[Mi-n]A3 
U2V2 



g 2 qi 

[n]qA2 



U 2 V! 



Ai 
U2V2 




M 

V2V1 









Ai 

U2V2 




A4 





Ai 

U 2 V 2 



9 2 Qi 

g 2 z2 




U3V1 



Q2Z2 





Aa 

U 2 V! 





(A.22) 



/ U2V2 
A0A4- 1 





A- l = 



g 2 giUlViV2 
g 2 A z 2 





[n] q U 2 V 2 ~g 2 \M x -n\ qi UlV x V2 
A0A2 1 g 2 A AiA^ 1 z : 





V 








Z-2 



U2V2 
A! 





AaA x 1 

_JJM o 

A>-4 4 

g 2 qiq2UlViVi [Nh-n] U 2 Vi 
-g 2 A Ai -AoA^ 1 









g2C/|ViV2i2 





Mgggj^fVWgg g 2 qiq2U'iV?V2 

-AoA^Ai -g 2 AoM 

[Nh-n\U 2 Vi 
A0A3 1 



-A A 2 







U2V1 

At 



q2UlVlV2Z2 

A 

g 2 giujv? , 

g 2 AoAi Mz 2 J 



(A.23) 



here we have defined &• 



gag 
9 

T 2 jr2 



and „4, 







-4.1 = bi^U^Vf - a 2 bi, 
A3 = a 2 ai 



[n] q A 1 A 2 + [Mi 



n 



, A3A4 where 



A 2 
Aa 



c 2 ciU% - ciVfc 2 , 
d\c 2 Vi - c 2 diU 2 . 



(A.24) 



B 



-[Mi-fci],g 2 B3 
[M 2 -k 2 ] q q 3 B7 
qiUM 




[ki] q q:jBj 
V1V2 



v\ 1 



U1V1 







gif/iVi 
-[Mi-fei],g 2 B3 









fc2] g <?3S7 



V 







[fci] 9 g36 4 
ViV 2 




[A/2— r 

qiUiVi 
q S B 5 
V1V2 




[k 2 ] q Bs 









_ 53S5 
B&_ 





<72#2 


[A/i-A;i] g 92B3 



9361 \ 



UxVi 






qiUiVi 


[M2-fc 2 ] g g3g7 

giC/iVi 




9365 



ViV 2 



(A.25) 



and we are using the shorthand notation q\ 

a 2 h, 



B 7 



b 2 a 2 U?V? 



a 2 a\ — a\a 2 

-2 



a 2 a 2 (l 



Mi 
2 



, Q2 

B 2 = ha 2 - 
Ba = c 2 c\V 2 
Be--*-^' 2 



2 , 



93 = g 



02&1, 

- ClC 2 V! 2 , 
(I2C2V1 - C2CI2U1, 

c 2 c 2 (C/i 2 - V?). 



fe 2 - 



M2 



and 



(A.26) 
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The matrix S^f 1 '^ 2 is defined as 

= I "\ -JLm ) > (A.27) 




where only first three rows of both and ' are taken. 

Rational limit. In the rational limit q — >• 1 the coefficients (A. 24) and (A. 26) acquire 
quite compact expressions 

Ai _ /"y r~g~{ x i ~ x i) (1-^2)72 



giving 



a a V Mi y M2 x x £ 2 71 

— = aA 2 = i x —\— — , A.28 

a Y M 1 V ^ 2 x 2 



and also 



a a Bp, = i- 

""l J '2 x 2 



x l x 2 x 2 



£2 . 
— - = a a 05 = z 



' 9 


I 9 Oi 


- xf) (l - x l xj) 72 


M\ \ 


I M 2 


^2 7i 



^3 2 ~„ . rv r~9~ ( x t - x t) ji^2 

— = a aft - 



5 \ Mi\ M 2 xjxj 

&7 2~ n -9 ( X l ~ X t) 72 / a on\ 

— = a z az3 8 = z— ^-p . (A.30) 

a M 2 x 1 x^ 

B Elements of the special cases of the S-matrix 
B.l Elements of the fundamental S-matrix 

The fundamental S-matrix for the space III acquires the following form, 

s|0V> = (J^^V) + (2o 1,0 )il^V> + (^ 1 1 '°)^ 1 v 2 ) + (^^V), 

§ |0V> = (^ 1 °' 1 ) 1 |0 1 ^ 2 ) + (J^ 1 )^ 2 /) + (jf 1 )?^ 2 ) + (^o 0,1 )?I^V), 

§|VV> = (^ 1 1 '°)||0V 2 > + (^o 1,0 )5^V) + (5r 1 1,0 )il^V 2 > + (^ 1,0 )^V), 

§ |^V> = (jf' 1 )*^ 2 ) + (iT 0,1 )^V> + (5f '^ilVV 2 ) + (^^V). (B.l) 
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In order to find these coefficients 2? it is sufficient to consider the first relation of (4.67) 
and its affine counterpart only. In fact, the constraints read as follows, 



'{G°p)1 {G°P)f 

{G°P)\ (G°P)l 

'{G°P)i {G°P)f 

(G°P)i (G°p)1 

'(G<*)\ (GJP)t 

(G°p){ (G op )l^ 

\g°p)\ {G°p)1" 

(G°p){ (G°p)1 



1(1,0) 



1(1,0) 



(0,1) 



1(0,1) 



1,0\ 



3f^)\ (jf' 1 )* 



oyl,0\l ( or 1 




itself^ ' 1 ) 1 , 



(G)l (Gf)2 

(G)i (5)2 

(G)J 

(G)? (G)2 
(G)f (G)§ 

'(G)! (G)g 



(1,0), 



(0,1), 



(1,0), 



(G)i (G)J 



(0,1). (B.2) 



It is easy to solve these relations for 2f and we find them to agree with [10]. 
completeness, we have listed the relations of our elements 3? to those of [10] 5 



For the 



(5i ' u )f (J^ 




2f^)\ {2°> x )\ 

0,1n5 



0,^5 
6, 



j^i' 1 )? (ir; 

<5V=0,l\i /dj«0,l^i 



^0,1^6 
•^0 



0,^6 
6, 




1 
1 

^12 



-F12 ' 

q+q- 1 

Di2-E 12 



q+q * 9+9 
q- 1 A 12 +qB 12 



q+q 



-i 



qF 12 



qC 12 

q+q- 1 

' qA 12 +q- 1 B 12 

q+q- 1 
_qZl£ 12 
q+q- 1 



q+q- 1 

q Di 2 +qE 12 

q+q- 1 

q~ 1 F 12 

q+q- x _ 

qDi 2 +q 1 E 12 

q+q- 1 



(B.3) 



B.2 Elements of the S-matrix Sqi 

Here we list the explicit forms of the coefficients of the matrix Sq\. 

Subspace II. First we give the coefficients of the matrix *3f in the case of a bound 
state scattering with a fundamental particle. There are four different combinations of the 
parameters k±,k2,n that contribute. Thus we have to consider the case where &2 = and 
k\ = n = k leading to 



(*2S 



fc,0\l 

1 



,h+k 



U2V2- 



"1 - x 2 zu - q 



Q-2k-l 



X 



212 - q Q ~ 



1=2 \Q - k] q x 2 - x% U2V2 71 



U{Vi 72 



l-Q _ 

' 2 aU 2 V2[x 1 



X-l Xn 



TV 



q- Q [k] q xt 



x\ )(x 1 x 2 - 1)7172 
- x 2 x± 7172 



^/\Q} q (x7-x+)(l 



a 



k,0^2 
2 



(9£S)? 



(0£l°l)2 



1 



•TV 



g 2 rJiVi ^1 



1 



"l -TV7 72 



Xo 71 



0. 



(B.4) 



2 x 



■ 'V ■ I V » .'V'* _ 



1 n t x 



Xi Xn 



5 We remind that our as parameterization is based on the one of [35] which are related to those of [10] by 
351 = SS — 1 ( a; [io] + 0- This point must be taken into account when performing the concrete comparison. 
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Next we have three elements corresponding to k 2 = 1 an d k\ + 1 = n = k giving 
x i ~ x 2 (<? 2< ' fc+1 - ) — q 2 ®)z\2 



(f£S)i 



Q u 2 v 2 - 



(Of/kfisi _ [Q - k - l]q 



zn - q 



1-1 



(»£S)? 



7i72 



■jQ-k-l 



V\Q\q( x i ~ x 2 )(!- x 1 x 2 )xj a 



qi+k Q Xl — xf 72 
\/\QU x i ~ x t^l 
(B.5) 



Then we have another three scattering entries for k 2 = and fci = n + 1 = k contributing 



AOM _ A+Qtt ^ x 1 ~ x 2 1 ~Q 



(»£S)! = ^Wa^ 



-2A; 



37 n 



?2l2 



1 + Q-2A; 



[A;], x 2 - U 2 V 2 7i 



xt UiVi 72 



(s£i°i)i 



fc,0\l 



(B.6) 



Finally, there is one element with ^2 = 1 and k\ = n = k — 1 providing the last element 



(B.7) 



- x 2 - g^i2 

Subspace III. There are 36 elements of the matrix 3f that need be determined. As 
mentioned in Section 5, it follows that (5.16) becomes 

S|A; — 1, l)f = 9 (\k, 0)f l +q\k-l, - qS\k, 0)f l . (B.8) 

Acting with the S-matrix on both sides of the equations (5.17) and using its invariance 
property allows us to express the elements of the S-matrix of the left hand side to the ones 
on the right hand side. Explicitly we find 



fc,0\lr 



-Q-2_ q) + ^kfi^i\Q fc + 



fe,0 



l [k-l] q [Q-k + 2] q tf 



[Q-2k + l\ 



k,0 \1 



[k]c 



'i j 



(^-1.1)3 =(5 j 



fc,0\3 



[k - l] q [Q - k] q q 2k ~ Q ~ 2 - q[k] q [Q - k + 1], 



+ (^£1)1 



[k]c 



[Q - kl 



( JT^ - ^'^)f = ( ^ K ' v 



, i0 , 3 [* -2] q [Q-k+ l] q q 2k ~Q- 4 + # - l] q [k -Q- 2] q + [2k - Q - l] q 

h — 



;[*-!]< 



(X 



fc-l,l\5 

h 



, ^5 ([^-l] g g 2fc ~ 3 ~ Q " - k + l} q 

'-y-^k h 



k,0 n6 



( j2° ) =(j%" ' 



[*] s 

([fc-l] gg 2fc - 3 -< 

[*], 



#] g )[Q-A: + l], 



(•Sfc 1 ' 1 )! =(Cl U l)3 



fe,0 si 



[fcyQ-fc + i] 9g 2fc -Q- 2 

[* - 1], 

" [ V^ 2]g + [ fc -g +i ] 



g[Q - fc] g 



. /^fc,0\l [Q~ k+ l] q 
+ l-^fc_iJ3- 



[fc - l] c 



-fc+n 



fc,0\3r 



2fc-Q-2 



+ (^l)3 



[k - 1], 



+ (^T' U )3 



MM Me 



[k-l], 
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h-^k-lhy q Q-2k+4 [k _ 1]q • [k _ l]q 
(^-1,1)5 =(a *0 ) 5 ([ g _ fc + l]qq 2k-Q-3 _ g[Q _ > 
(^-1,1)6 = ( a *0 ) a ([ g _ fc + 1]gg 2 fe -Q-3 _ ?[Q _ k]) _ 



fc,0\3 



fc-2] ff [Q- fc + 1], gtfcyfe-Q + llg-^pfc-Q-l], 



+ 



+ (3f'°)§ , 
(B.9) 



Finally, the remaining elements are 
a ( x i ~ x t)( x 2 — 4) 



fc,0\l 
fc 



«L°i)I 



(3i 



fc,0\5 
fc,0 x6 



( ^, 0) 1 



[^]<?(4 4 



Q 
• 2 



UlVl 7i72[£W[QKl " x i x 2~)( x i - 4)^ 

q fo[ - 4)(x 2 - 4)[g(£ + x^)(gx+ + 1) - (g + x+)(£x+ + 1)] 



7i72(C 2 - 1)V [Q](l - Si x 2 )Oi - 4)<? 



71 [g ~ k] q (aj - x+)[g(g + x^)(g + gj) - (gx^ + l)(gx+ + 1)] 
" 72 <P vtt (C 2 -l)(l-xrx2)C7 1 y 1 (xr-x+) 
.71 [* - 1], (xj - 4)fe(£ + xi )(£ + x 2 ) - (^ + l)(£x+ + 1)] 



72 



(^ 2 -i)(i-s^2)^i^r-^ 



. (x+ - 4)[(^r + i)(c4 + 1) - ?(£ + *r)(£ + ^ )] 



Q+l 
2 



(£ 2 - l)(xrx 2 " - l)(xi - x+)i7 1 F 1 C/ 2 ^- 

_z 12 (x2 - x+)(x+(£*2 + i)(ex+ + 1) - + x 2 )(e + 4)) ^ 



(e 2 - - i)K - 4)<?^ ' 

- x+)[q - fc] Q [(exr + i)(£c+ + 1) - y 2 2 (^ + xr)(e + 4)] , 



-2 ^1 



(£ 2 - 1)4^12 [Q]<?(4 x 2 -l)(zi -4)? Q 
xj~(x 2 x]*~ — 1)(4 ~~ X2)q k ~ 2 ® 

x f( x l x 2 ~ ^)( x l ~ 4) 

x^[k} q q- k (x^ - x+)[(gx7 + l)(£c+ + 1) - F 2 2 (g + x7)(g + 4)1 

(1 - £ 2 )x^l 2 [Q] ? (l - X5"X2")(X5" - X^) 

, Tf arWgg"^ 1 ^ + ~ fc Ug4(£ + ^2 ) ~ x 2 (£ + 4)] 
"(£ 2 - l)xfx$[Q] q (l - x\x^ )(xj" - x£) 

_ 7l 2 xr[A: - l],[A;]g- fe - 1 (^+ + l)[gx+(£ + x 2 ) - x 2 (£ + x+)] 



(£ 2 - l)x+x+[Q] 9 (l - x x x 2 )(x x - xj) 



: 7i72(^ 2 2 - i)gr^2 [ fc ]<?g~ Q ~ 5 (4 - 4)(£ + 4) (£4 + 1) 
a(£ 2 - i)44v / [Ql?( 1 ~ x r x 2 X^r - 4X^2 - 4) C/ 2V 2 ' 

7i72^r [/c] g(? ~ Q ~2(a;+ - x£) 



"4 VTOl? ( x 1 x 2 - !)( x l -4)^2^2 

qg(l - y 2 2 )(xr - x+) 2 (g + x 2 )(gx 2 - + 1) 

7i 2 (£ 2 -i)[Q],(xrx 2 --i)(xr-x+) ' 

q(l - F 2 2 )g fc ~ 3 (xr - x+) 2 (g + x 2 -)(gx 2 + 1) 

7 2 (£ 2 -i)[QUi-xrx 2 )(xr-4) 
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/j»fc,cK3 = q k ( x i -x 2 )(x x x+ - l) | 

( x l x 2 ~ l)( x l ~ x 2 ) 

x 2 -M g (xr - x twKi x i + i)M + 1) - (j + ^r)(c + 4)} 

(£ 2 - l)x^i 2 [Q] 9 (x^x 2 - - xj) 

_ x 2 [fe - i] g r fc (^r - + i)M tmil x D(e + ^)] 

(£ 2 - l)x^l 2 [Q]g(x5"X2 - - X%) 

, ^k,o-s5 _ I Q ( v 2 ~ l)(a?i - a?i")(l - a?i «2 )(£ + x i)^ x 2 + 



[Q]g (C 2 - x 2 -l)(x 1 -sj)(z 2 -^2) 
C Yang-Baxter equation 

In this section we briefly summarize some details on the checks of the Yang-Baxter equation 
(YBE) that we have preformed for the bound state S-matrix. 

Let us first focus on subspace I. The block X governing the scattering in this subspace 
is required to satisfy YBE on its own right. Thus we need to consider the following 
scattering sequences, 

|0, 1, fei, ki) <g> |0, 1, k 2 ,k 2 ) <S> |0, 1, k 3 ,k 3 ) > |0, 1, mi, mi) <g> |0, l,m 2 ,m 2 ) <8> |0, 1, m 3 ,m 3 ) 

(C.l) 

which give the explicit form of the YBE in subspace I, 

n=0 

fc2+fe3 

£ ^rr n ' n (^^2)^^ +fc3 - n (^i,.3)^ 2 ' fc3 (^^3). (c.2) 

n=0 

We did not attempt to prove this identity in full generality, but we did check it for a large 
set of different values of the parameters ki,rrii and bound state numbers and found it to 
be perfectly satisfied. 

In a similar way, this approach for checking YBE may be extended to include the 
subspaces II and III. For example, acting with Si2S± 3 S2 3 — S23S13S12 on the states of the 
following form, 

|0,0,fci,Ji) ® |0,0,fc a ,Z a ) ® \0,l,k 3 ,k 3 ), (C.3) 

will result in a plethora of different types of states and the coefficients will depend on all 
three scattering blocks X ', *3f , 3?. Due to the large size of these expressions we have not 
spelled them out explicitly here. Nevertheless we have explicitly computed, for different 
values of the parameters, various matrix elements of the YBE that include states from all 
three subspaces that in general may be written as 

( out-state I YBE | in-state ) . (C.4) 
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We have performed the checks for a wide range of numerical values of the representation 
parameters and in each case it proved to be compatible with the YBE. 
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